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The statistical mechanics of particles with shapes on a one-dimensional lattice is investigated 
in the context of the s = 1 Ising chain with uniform nearest-neighbor coupling, quadratic single- 
site potential, and magnetic field, which supports four distinct ground states: | titi "')> I 00 '''}) 
ltt''')ilt t '''}- The complete spectrum is generated from each ground state by particles from 
a different set of six or seven species. Particles and elements of pseudo-vacuum are characterized 
by motifs (patterns of several consecutive site variables). Particles are floating objects that can be 
placed into open slots on the lattice. Open slots are recognized as permissible links between motifs. 
The energy of a particle varies between species but is independent of where it is placed. Placement of 
one particle changes the open-slot configuration for particles of all species. This statistical interaction 
is encoded in a generalized Pauli principle, from which the multiplicity of states for a given particle 
combination is determined and used for the exact statistical mechanical analysis. Particles from all 
species belong to one of four categories: compacts, hosts, tags, or hybrids. Compacts and hosts find 
open slots in segments of pseudo-vacuum. Tags find open slots inside hosts. Hybrids are tags with 
hosting capability. In the taxonomy of particles proposed here, 'species' is indicative of structure 
and 'category' indicative of function. The hosting function splits the Pauli principle into exclusion 
and accommodation parts. Near phase boundaries, the state of the Ising chain at low temperature 
is akin to that of miscible or immiscible liquids with particles from one species acting as surfactant 
molecules. 



I. INTRODUCTION 

The constituent particles of condensed matter are 
subject to strong interaction forces. Inspired by the 
search for separable coordinates in Hamiltonian mechan- 
ics, much progress in the understanding of material prop- 
erties has resulted from finding collective modes that are 
much more weakly interacting. The (quasi-)particles as- 
sociated with these modes are generated out of vacua rep- 
resenting particular solid states. Exotic exclusion statis- 
tics are one manifestation of this diversity. 

In the study of magnetic insulators, commonly mod- 
elled as lattice systems of interacting spins, a variety of 
particles have been identified to describe their statistical 
mechanical properties under specific circumstances [T]- 
3 . In each case the physical vacuum (ground state) is 
configured as the pseudo-vacuum for particles from one 
or several species representing collective excitations. In 
general, these particles scatter off each other elastically or 
inelastically. Particles may decay or form bound states. 
An exact solution is almost always out of reach. 

To the rare exceptions belong the completely integrable 
spin chain models. In these systems, the particles only 
interact via two-body elastic scattering. The kinematic 
constraints of one dimension combined with the dynamic 
constraints of a factorizing S'-matrix impose conservation 
laws that guarantee an exact solution (via Bethe ansatz, 
for example) [4]-[6]. Many integrable spin model systems 
have one or several Ising limits (of commuting spin oper- 
ators) , where some of the key particle properties such as 
their exclusion statistics are left intact. 

Here we propose a taxonomy of free particles in Ising 
chains. We classify particles according to both struc- 
ture and function. We associate structures with particle 



species and functions with particle categories. It appears 
that Ising chains produce particles of a seemingly unlim- 
ited diversity of structures but of only four categories. 

It remains to be seen how robust some of the statisti- 
cal mechanical properties of the particles identified here 
are when non-commuting terms are added to the Ising 
Hamiltonian. The insights reported in this study are not 
limited to magnetism. Ising systems are applicable to 
phenomena in other areas of physics and beyond. The 
methodology introduced here is adaptable to applications 
beyond Ising models. 

We consider the spin-1 Ising chain with periodic 
boundary conditions, 



H 



i=i 



J Si S i 



I J l + 1 



D{St 



hSf 



Sf = 0,±1, (1) 



where the J-term denotes a nearest-neighbor coupling, 
the D-term a quadratic potential, and the ft.-term the 
Zeeman energy. In higher dimensions this Hamiltonian 
coincides with the thoroughly studied Blume-Emery- 
Griffiths model [7j. The statistical mechanics of the spin 
chain per se is of limited interest and most conve- 
niently worked out via the transfer matrix method [5]. 
We employ this textbook model as a vehicle to develop 
an idea of much wider scope. 

The spectrum of H is expressible in the form of prod- 
uct states, |o"i<72 • • • crjv) with oi =t, The zero- 
temperature phase diagram features six phases as shown 
in Fig. [l}a) [5] : one antiferromagnetic phase <I> a > one sin- 
glet phase $s, two ferromagnetic phases $f±, and two 
plateau phases &p± [10]. The physical vacuum (ground 
state) is unique in &f±, &s & n d twofold degenerate in 
<£ A , $ P± . 
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FIG. 1: (a) Six T = phases $a (antiferromagnetic), $f± 
(ferromagnetic), $s (singlet), and <J>p± (plateau) for J > 0; 
(b) three phases at ft = 0. The text refers to regions in 
parameter space by the same names. 



In Sec. [TT] we configure each physical vacuum as the 
pseudo-vacuum of a set of dynamically free but statis- 
tically interacting particles. The statistical mechanical 
analysis, carried out in Sec. |III| yields the exact popu- 
lation densities of particles. In Sec. |IV| we interpret the 
low-T state near phase boundaries as that of miscible or 
immiscible liquids. Section [V] further elaborates on the 
taxonomy of particles. 



II. COMBINATORIAL ANALYSIS 

It will suffice to consider four sets of particles owing to 
spin-flip symmetry. The two sets generated from twofold 
pscudo- vacua (§a,®p+) have shared attributes. The 
two sets generated from unique pseudo-vacua (<1>5, $f+) 
share different attributes. Each choice of pseudo- vacuum 
can be used in the entire parameter space, but only in 
one region does it coincide with the physical vacuum. 

Particle species and elements of pseudo-vacuum are 
characterized by motifs comprising two, three or four con- 
secutive site variables ct/. Motifs interlink in one shared 
site variable. Particles are subject to mutual statistical 
interactions as captured in the generalized Pauli principle 
proposed by Haldane [TT]. The relation, 
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expresses how the number of slots available to one parti- 
cle of species m is affected when particles of any species 
to' from a given set are added. This idea has proven 
to be surprisingly versatile with a wide range of applica- 
tions and a high degree of flexibility regarding the notion 
of particle species [31 IT2TH4] . 

In a chain of TV sites the number of many-body states 
containing a specific combination of particles from all 



species is determined by a multiplicity expression inferred 
from ([2]) via combinatorial analysis [UJ . The exact struc- 
ture of the multiplicity expression varies between appli- 
cations. We have developed a template that can be em- 
ployed for all situations considered here: 
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where d m counts the number of open slots for particles 
of species m in the presence of N m i particles from any 
species to'. The permissible combinations {N m } produce 
integer- valued d m . 

Each application is characterized by specific values 
of M (number of particle species), n pv (multiplicity of 
pseudo- vacuum), A m (capacity constants), a m (size con- 
stants), and g mm i (statistical interaction coefficients). 
These specifications will be tabulated for each situation 



under scrutiny. The prefactor in ( 3a I is specific to the use 
of periodic boundary conditions. It is a generalization of 
a modification factor worked out by Polychronakos [T3] , 
The energy of a many-body state only depends on its 
particle content, not on the particle configuration: 
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where E pv is the energy of the pseudo-vacuum and e m 
the particle energy (species to) relative to the pseudo- 
vacuum. There is no interaction energy. Each phase 
in Fig. [TJa) is the pseudo-vacuum of one set of particle 
species or a solid of one species of particles (with e m < 0) 
from the same set. 

All particle motifs occupy at least one bond on the 
lattice. Adding particles of any species, therefore, dimin- 
ishes the number of open slots for more particles. Nor- 
mally we have g mm ' > 0, which enforces Ad m < in ^ 
whenever a particle is added. The limited space for parti- 



cles on the finite lattice is signalled in ( 3a I by a vanishing 
binomial factor caused by the factorial of a negative d m 
in the denominator. 

However, we shall find some g m m' to be zero or nega- 
tive. A vanishing g mm i means that placing a particle to' 
has no direct effect on the options for placing a particle 
to thereafter. A negative g m m' turns the exclusion prin- 
ciple into an accommodation principle. Placing a particle 
to' increases the options for placing a particle to there- 
after. In these instances, space limitations on the lattice 
are enforced by indirect effects on the d m . The statisti- 
cal mechanics of a macroscopic system only depends on 
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as will be demonstrated 
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A. Antiferromagnetic pseudo-vacuum 

In region the ground state is the twofold degener- 
ate Neel state, | tltl ■■■), I Itlt ■■■), implying n pv = 2. 
We have identified a set of M = 6 species of dynamically 
free particles excited from this state [16]. The motifs 
and statistical interaction specifications are summarized 
in Tables |T] and [TT] The full spectrum of product states 
\<J\ - ■ ■ on) is equivalent to all possible particle configura- 
tions. The number of states with given particle content 
{N m } and energy Q is determined by (3a I. 



TABLE I: Specifications of M — 6 species of particles excited 
from the Neel state (n pv = 2), | tltl ' ' ■ ). I It It ■ ■ ■ ) : motif, 
category, species, energy (relative to pseudo- vacuum) , spin, 
capacity constants, and size constants. Segments of £ vacuum 
bonds tli It have energy £(D — J). At h 7^ the entries of 
e m must be amended by —s m h. 
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Particles m 7^ 3 and vacuum bonds are placed side by 
side onto the lattice. They interlink in one shared site 
variable as in these examples: 
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(5) 



Particle 3 can only be placed inside particles 4,5,6 as in 
the following examples: 

toot =4 + 3, t°°°l =6 + 3 + 3, 

t o t 00 It 00 t = 6 + 5 + 3 + vac + 4 + 3. (6) 

The six particle species naturally divide into three cat- 
egories: compacts, hosts, tags. Particles m = 4, 5,6 are 
hosts to tag m = 3. The compacts m = 1, 2 do not host 
nor are they being hosted. 

The motifs of compacts 1,2 are to be interpreted as 
consisting of one bond and half of the site at either end. 
In the hosts 4,5,6 we count two bonds and two sites (the 
site in the middle and half of each outer site) . Tag 3 has 
one bond and one site. The site on the right is counted 
fully. The site on the left is counted as interior site of 
the surrounding host or as right-hand site of another tag 
to its left. With these rules each particle has a definite 
energy e m and a definite spin s m . 



TABLE II: Statistical interaction coefficients of the particles 
excited from the Neel state (or plateau state.) 



Qm m ' 


1 


2 


3 


4 


5 


6 


1 


1 

2 


1 

2 


1 

2 





1 


1 

2 


2 


1 

2 


1 

2 


1 

2 


1 





1 

2 


3 











-1 


-1 


-1 


4 


1 
2 


1 

2 


1 
2 


1 


1 


1 
2 


5 


1 
2 


1 

2 


1 
2 


1 


1 


1 
2 


6 


1 


1 


1 


2 


2 


2 



Vanishing and negative g mm ' are indeed realized in Ta- 
ble [n] The placement of one host, m! — 4, does not 
affect the number of open slots for placing a compact, 
m = 1; hence 514 = 0. Likewise, the options for placing 
a tag, m = 3, are unaffected by first placing a compact 
m! = 1, 2, or another tag; hence = g^2 = .933 = 0. 

If a host m! = 4, 5, 6 is placed into any open slot then 
the options for placing a tag m = 3 increases by one; 
hence 334 = 935 = .936 = — 1. Tags have vanishing capac- 
ity constants: A 3 = 0. Open slots for tags are created by 
hosts. 

The two distinct motifs of host 6 do not compromise 
the indistinguishability of particles belonging to the same 
species. At most one of the two motifs can be placed 
into any particular slot. Distinguishable motifs of indis- 
tinguishable particles cannot exchange positions. It is 
possible to split species 6 into two species with one mo- 
tif each. This introduces a mutual statistical interaction 
between the two new species and modifies the statistical 
interactions with all other species (TTJ [IS] . 

In region <S>a all e m are positive. Here the physical 
vacuum is the pseudo-vacuum. Elsewhere in parameter 
space we have e m < for at least one species. There the 
physical vacuum is a solid formed by the particles with 
the lowest (negative) energy density. We can reconfigure 
any one of these solids as the pseudo-vacuum for a dif- 
ferent set of particles. In Ref. [16] these particles were 
interpreted as soliton pairs in generalization of the well- 
known classification used in the context of spin-^ XX Z 
chains |19) . 



B. Singlet pseudo-vacuum 

In region $s the ground state is unique, | o o • • - o), 
implying n pv = 1. It is the pseudo- vacuum for a set of 
M = 7 particle species. The specifications are listed in 
and IV The multiplicity expression ^ again 



Tables III 



accounts for the full spectrum. 

There are no compacts among the seven species. Parti- 
cles 1,2,3,4 are hosts to particles 5,6,7. The latter behave 
like tags in some respects. However, their ability to co- 
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TABLE III: Specifications of M = 7 particles excited from 
the singlet state (n pv — 1), | o o ■ ■ ■ }: motif, category, species, 
energy (relative to pseudo- vacuum) , spin, capacity constants, 
and size constants. Segments of vacuum bonds oo have zero 
energy. At ft 7^ the entries of e m must be amended by —s m h. 
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TABLE IV: Statistical interaction coefficients of the particles 
excited from the singlet state (or spin-polarized state.) 
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exist inside the same host implies a hosting capability 
among them. We have named them hybrids. Particle 7 
has two motifs, which again causes no complications. 

Each host is an island in the singlet pseudo-vacuum. 
The size of an island depends on its hybrid content. 
When a hybrid is placed inside a host the first site vari- 
able (on the left) of its motif is shared with an interior 
site variable of the host or any matching site variable 
of a hybrid already present inside that host as in these 
examples: 



o ft o = 1 + 5, o f4-t4- o = 3 + 7, 

{°ntt°,°mt°} = 1 + 5 + 7. 



(7) 



In the last instance, two states correspond to the same 
combinations of particles, which is accounted for by the 
multiplicity expression (|3]). 

In the course of adding a hybrid, one bond of the host 
or of a hybrid already present is broken up and replaced 
by an identical bond to the right of the newly placed 
hybrid. Each host has at most one open slot for any par- 
ticular hybrid. The second example in ^ is no exception 
because it makes no difference whether one or the other 
motif of hybrid 7 is placed inside host 3. 

Whereas a given combination of particles may be 
placed in multiple ways on the lattice, each product state 



has a unique particle content as expressed by the list 
{N m } in (3a) and Q. The rules that ensure that each 
bond and each site are counted exactly once are the same 
as in Sec. Ill Al 

The g mm ' listed in Table |TV| again include some vanish- 
ing or negative coefficients. Not all occurrences have the 
same interpretation. Hosts ml = 3, 4 can accommodate 
all hybrids m = 5, 6, 7 whereas hybrid m = 5 (6) can- 
not be accommodated by host m! — 2 (1) without the 
presence of at least one hybrid m — 7. This explains the 
non-positive coefficients g mm ' between hosts and hybrids. 

Ambiguities arise when we aim to determine the g mm i 
among hybrids. Hybrids m = 5, 6 alone do not share 
open slots for placement, implying 1756 = g§§ = 0. New 
situations present themselves when hybrids m — 7 are 
added to the mix. All hosts are capable of accommo- 
dating hybrids m = 5, 7 or m — 6,7 or m — 5, 6, 7 in 
multiple combinations. 

Consider a single host, say m = 1, with several hybrids 
embedded. If all hybrids are of the same species, e.g. 
m! = 5 or m! = 7, then the multiplicity is dictated by the 
size of the expanded host. This information is encoded 
in the g mm > between host m and hybrid m! . 

However, if host m = 1 accommodates hybrids m! — 
5, 7 in several different sequences, then that contribution 
to the multiplicity is encoded in the statistical interac- 
tion coefficients between the two hybrids. Since there is 
no hierarchy between the two hybrids, the encoding of 
that information can be implemented in #75 or in g^. In 
Table [TV] we use 557 = g&j = —1, 575 = gi§ — 0. Different 
choices produc e the same multiplicity. The hybrids (like 
the tag in Sec. II A) have zero capacity constants A m . 



C. Ferromagnetic pseudo-vacuum 

Here we pick the spin-polarized state, | tt • • • ), as the 
pseudo-vacuum for a set of particle species that again 
generate the full spectrum. This state coincides with the 
physical vacuum in region $ p + . The specifications for a 
set of M — 7 particles can be obtained, in this instance, 
by reasoning as follows. With the transcription 



{°,t4}s = {t, + 4+F+ 



(8) 



between site variables we produce a one-on-one map- 
ping between Ising product states that maps the singlet 
state I o o • ■ • ) onto the spin-polarized state | ft " " " ) • 
Given that the singlet state is the pseudo- vacuum of the 
particles listed in Table |III| the transcription yields a 
corresponding set of particles generated from the spin- 
polarized state. These particles are listed in Tabic |V) 

Hosts map onto hosts and hybrids onto hybrids, leav- 
ing the A m and a m invariant. The g mm ' keep the values 



listed in Table IV The mapping changes the specifica- 
tions e m , s m . Given the nonzero spin-polarization of the 
pseudo-vacuum, the quantum number s m is not a spin 
in the usual sense, but its use enables us to write the 
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TABLE V: Specifications of M = 7 particles excited from 
the spin-polarized state (n pv = 1), j ft ' ' ' ) : motif, category, 
species, energy (relative to pseudo- vacuum) , 'spin', capacity 
constants, and size constants. Segments of £ vacuum bonds 
tt have energy £(D + J). At h ^ the entries of e m must be 
amended by —s m h. 
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magnetic-field dependence of the particle energy in the 
Zeeman form, —s m h. 

D. Plateau pseudo- vacuum 

The partially magnetized plateau state, | f t ' ' ' )> 
1° t t ' ' ' )) which is the physical vacuum in region $p+, 
is selected here as the pseudo-vacuum for the fourth set 
of particle species. These particles can be found with the 
transcription 

MiU = {U,°}p+ (9) 

between site variables, which produces a one-on-one map- 
ping between the Neel state and the plateau state. The 
specifications of the M = 6 particles transcribed from 
Table |T] are listed in Table |VI| Only the motifs and the 
entries of e m , s m are different. The g mm i remain as stated 
in Table [TT] The transcription does not affect the cate- 
gories. 

TABLE VI: Specifications of M = 6 particles excited from 
the plateau state [n pv = 2), ]^o-^o--.}, |ot~o"t---): motif, 
category, species, energy (relative to pseudo- vacuum) , 'spin', 
capacity constants, and size constants. Segments of I vacuum 
bonds t °, ° t have energy ID/2. At h 7^ the entries for e m 
must be amended by — s m h. 
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The structural correspondence between particle species 
as determined by transcriptions (JsJ) and ^ does not im- 
ply a functional correspondence because some functions 
depend on the relative energies as will be discussed in 
Sees. |IV] and El 



E. Species and categories 

The following comments are based on wider evidence 
than reported in this work, including evidence from Ising 
chains with spin s = | , | and with next- nearest-neighbor 
coupling, where sets of particles with up to 17 species and 
motifs of up to six sites have been found [20 . 

(i) The template ^ for the multiplicity expression 
with specs np V ,A m ,a m ,g mm i characterizing the pseudo- 
vacuum and all particle species excited from it has proven 
sufficiently robust to hold for all cases investigated. 

(ii) The classification of particles into the four cate- 
gories of compacts, hosts, tags, and hybrids has proven 
useful and comprehensive within the range of our explo- 
rations. Compacts and hosts find open slots on segments 
of pseudo- vacuum, tags and hybrids inside hosts, hybrids 
with hosting capability for different hybrids. 

(iii) The capacity constants A m are extensive (oc N) for 
compacts and hosts. The value of A m reflects the holding 
capacity of the pseudo- vacuum for particles of these two 
categories. It depends on the size of the motif(s) of a 
given species and how closely particles of that species can 
be packed. The pseudo- vacuum has no holding capacity 
for tags and hybrids, implying A m = 0. Holding capacity 
for tags and hybrids must be generated by hosts. 

(iv) The statistical interaction coefficient g mm i reflects 
the direct impact on the open slots available to parti- 
cles of species m if one particle of species m' is added. 
Negative coefficients g m m' represent an accommodation 
principle. The presence of hosts m! opens up slots for 
tags or hybrids m. However, for every g mm > < there 
is a g m ' m > 0, which, in combination, limit the holding 
capacity for hosts and tags (or hybrids). 

(v) The limited room for particles of any species is 
signalled in the multiplicity expression ([3| by a vanishing 
binomial factor as caused by factorial of a negative d m 
(interpreted as a T-function) in the denominator. 

III. STATISTICAL MECHANICAL ANALYSIS 

The statistical mechanics of the four sets of particles 
identified in Sec.[H]is amenable to a rigorous analysis that 
starts from the multiplicity expression ([3|. The method 
was developed by Wu [21] for a generic situation (see also 
Refs. [221 E3|). 

The range of applications is surprisingly wide. It in- 
cludes an alternative to the thermodynamic Bcthe ansatz 
|14j , applications to particles in real space and in recipro- 
cal space [21 [HI [T31 [2H [5S] , and applications to systems 
in higher dimensions [26] . Each application requires a 
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certain degree of adaptation of Wu's generic result. The 
application to interlinking particles with shapes is no ex- 
ception. 

For given energies e m and statistical interaction spec- 
ifications A m , g mm i , the grand partition function of any 
one of the four sets of particles has the form 



z = n ( +Wn 



m— 1 



(10) 



where the (real, positive) w m are the solutions of the 
coupled nonlinear algebraic equations, 
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With the w rn thus determined, the average number of 
particles, (N m ), are the solutions of the coupled linear 
equations [27] . 
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Note the sensitivity of Eqs. ( 10 ) and ( 12 ) to parti- 
cle category. Tags and hybrids, which have A m = 0, 
do not contribute their own factor to Z. Their contri- 
butions manifest themselves indirectly via their hosts. 



Equations (12 1 with indices m pertaining to tags or hy- 



brids are homogeneous, which makes the tag and hybrid 
populations conditional on host populations. 

The solution of Eqs. ( 10 )-( 12 ) covers a wider territory 
than the transfer-matrix solution of the Ising chain (I]) 
because we are free to assign any values to the particles 
energies e m [55] • We choose the tabulated values per- 
taining to the three-parameter Ising model for mere con- 
venience. The ordering tendencies and phase boundaries 
are then as shown in Fig. [T] 

In the following we determine the average population 
densities, (N m ) = (N m )/N, of particles from all species. 
These quantities inform us about the roles of particle 
species in the ordering tendencies and about the signif- 
icance of particle categories in the behavior near phase 
boundaries. 



A. View from $a 

The statistical mechanical analysis of the particles 
identified in Sec. |II A| begins with the solution of 
Eqs. (111. The excitation energies e m including Zeeman 
h, are taken from Table \L 



term, — s 

interaction coefficients from Table [IT 



and the statistical 
The six equations 



with these specifications can be simplified into 



= e K ' 



1 + W2 
1 + W 3 



»'G 
W 3 

WiW 3 _ K] W 2 l + W A 

1 + W4 ' W\ 1 + w$ 



_ e K., + K D +H 



1 + W2 _ 2H 

1 + Wi 6 ' 



= 1, 



1 + w 6 



= 1, (13) 



where Kj = J/k B T, K D = D/k B T, H = h/k B T are 
scaled parameters. Expression ( 10 1 for the grand parti- 
tion function becomes 



Z 



(1 + w 3 )e 



Kr 



N 



(14) 



Equations (13) reduce to the cubic equation for w 3 



The solution at h = is sufficiently compact to be cited: 
l + w 3 = c~ Kd cosh Kj + - 




e- K ° cosh Kj - - ) + 2e- R o . ( J ; 



We calculate the dependence of the population densi- 
ties (N m ) on the reduced parameters D/J and h/J at 
a fixed value of reduced temperature k B T/J by solving 
Eqs. (12) for given w m . 

In the high-T limit we have w\ = w-i — w 3 — wq = 2, 
W4 _ = W5 _= 3, and (JV X ) - (N 2 ) - (N 3 ) - (N 6 ) - §, 
(Ni) — (N5) — jg, implying Z — 3 N , the total multiplic- 



U/ = \Jv 5 ; = 
ity of states. 

Taking into account the sizes and multiplicity of motifs 
the population densities must satisfy the constraint, 



(Ni) + (N 2 ) + (N 3 )+2[(N 4 ) + (N 5 ) + (N 6 )j+ (N vac ) = 1, 

(16) 

where (-/V vac ) is the density of vacuum elements with mo- 
tifs Relation ([16]) holds at any T. For T = 00 we 
infer (N vac ) = §. 

As T is lowered toward zero, particles from all species 
with e m > are gradually frozen out and particles from 
some species with e m < are crowded out by parti- 
cles from the species with the lowest negative energy 
per bond. Inside any region in the parameter space of 
Fig. [lja) the state of the system in the limit T — > 
(physical vacuum) has either been emptied out of any 
particles (pseudo-vacuum) or has turned into a solid of 
one species of particles. The former is the case in region 
$^4, the latter in all other regions [50] , 

Compacts m = 1 have negative energies (ei < 0) in 
region Hence their high density in this region as is 

evident in Fig. [2](a). Phase is a solid of interlinked 
compacts 1. Tags 3 have negative energies in the region 
$s. That is reflected in Fig. [2jb) . These particles solidify 
at T = in that region. 

The plateau phase $p+ is a solid of interlinked hosts 
4 (with £4 < 0). The high density of these hosts in the 
narrow region at nonzero T is illustrated in Fig.[2]Jc). 
Hosts 6 have positive energy for any value oi h/J and 
D/J. Their average number remains small everywhere 
but reaches a flat maximum at h/J ~ D/J ~ as is 
illustrated in Fig. [2jd). 

The plots for (A3) and (Ne) are reflection symmetric 
with respect to the line h/ J = 0. Plots for ( N2) and (A5) 
would be mirror images of the plots for (Nx) and (N4), 
respectively. 
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FIG. 2: Average numbers (per site) (iV m ) of (a) compacts 
m = 1, (b) tags m = 3, (c) hosts m = 4, and (d) hosts m = 6 
versus D/J and h/ J at ksT/J = 1. The dashed lines indicate 
the phase boundaries as in Fig. [IJa) 



At the five phase boundaries that border region $a 
we have e m = for exactly one species m. In all but 
one case that species is a compact or a host. Particles 
from these categories interlink directly with elements of 
pseudo-vacuum and thus produce a state of high degen- 
eracy. 

The &a — &s phase boundary is the exception. Here 



the particles with e m = are tags. They do not interlink 
with elements of pseudo-vacuum. They exist only inside 
hosts but all hosts have e m > 0. The tags m — 3 by 
themselves can only form one zero-energy state, the sin- 
glet singlet state | o o • • • ) . The degeneracy remains low 
(threefold) . 

At the four phase boundaries not bordering region <$>a 
two species of particles have equal (negative) energy per 
bond. In two cases (3>s — $p±) the two species are a host 
and a tag, either m — 4, 3 or m = 5, 3. In the other two 
cases ($p± — 3>f±) the two species are a compact and a 
host, either m = 1, 4 or m = 2, 5. In all these cases the 
particle species involved interlink either side by side or 
one inside the other in many different configurations. 



B. View from $s 



The analysis for the particles identified in Sec. |II B| 
proceeds along the same lines. Equations (111 for this 
set of particles can be simplified into 



1 + w e 
l + w 5 ~" 
(1 + wi)w 6 
(1 + w 6 )w 7 

(l + w 2 )w 5 



e 



w 5 w 6 
1 + W7 

2Kj-K D -H 



-,K, 



1 + W4 



WlW 2 W 6 



(1+W2)(1+W3) 



(17) 



They reduce to a cubic equation for w$ (quadratic for 
h = 0). Expression (10 1 for the grand partition function 
reduces to ED 



Z = 



(1 + w 5 )e 



.Kj-K D +H 



(18) 



The dependence of the (N m ) on D/J and h/J at fixed 
ksT/J as inferred from Eqs. (12) are (selectively) de- 
picted in Fig. [3] 

At T = 00 we have the values w\ = |,tU2 = 1^3 = 4, 
w A = 3,w 5 = w e = 2, w 7 = 3; (Aq) = (N 2 ) = (N-) = ^ , 

(JV 3 ) = (A 4 ) = 2V, (^s) = (N 6 ) = I- At finite T the 
population densities are low for all particles in region 
Here all particles have positive energy, the ground state 
is free of particles, the physical vacuum is the pseudo- 
vacuum. In Sec. HI A| the same physical vacuum was a 
solid of tags from a different set of particles. 

Hybrids 7 have negative energy in region $^4, causing 
a high population density there as is evident in Fig. |3ja) . 
Phase $^4, which was the pseudo- vacuum of choice in 
Sec. |III A| now is a solid formed by these hybrids. 

In all other regions the ground state is a solid of one 
species of particles from both views. In regions $p± the 
particles with dominating population densities at low T 
are now hybrids 5 or 6. The contour plots (not shown) 
look very similar to Fig. 2ja) and its mirror image. 

Hosts 1 and 2 have a high population in regions &p± 
and solidify into plateau phases as T — > 0. The contour 




partition function for this case reads 



FIG. 3: Average numbers (per site) (Nm) of (a) hybrids m = 7 
and (b) hosts m — 3,4 versus D/J and h/J at fcsT/J = 1. 



plots for (Ni), (N2) look very similar to Fig. [2jc) and 
its mirror image. Hosts 3 and 4 tend to be frozen out or 
crowded out by other particle species in all regions. Their 
population densities remain low as shown in Fig. [3jb) . 

At two of the three phase boundaries that border re- 
gion $5 the particle species with e m = are hosts 
(to = 1,2), which interlink with elements of vacuum, and 
at the third phase boundary they are hybrids (m = 7), 
which do not. This explains the high degeneracy at the 
$ s — Qp± boundaries and the low degeneracy at the 
$5 — <&a boundary from the vantage point of the sin- 
glet pseudo-vacuum. 

At the remaining phase boundaries we again have two 
species of particles with equal and lowest (negative) en- 
ergy per bond. In all cases particles of the two species 
interlink and thus produce a high degeneracy albeit not 
for the same reason. In four cases they are hosts and 
hybrids. In two cases they are two species of hybrids. 



Views from <E>f± and $p± 



(1 + w e )e 



-2H 



N 



(20) 



The results for the (N m ) as inferred from Eqs. ( 12 ) con 



firm what is readily predicted by analogy to the previous 
cases. At low T we have a low density of any particle 
species in region $f+ because the physical vacuum of 
this region is the pseudo- vacuum. In regions $ p + , <&a 
we find high densities of hosts m = 1,2, respectively, 
and in regions &s,&f-,$p- high densities of hybrids 
to = 5,6,7, respectively. Hosts m — 3,4 do not reach 
significant densities in any region. 

At all phase boundaries except one the particles 
present in high density interlink directly among them- 
selves or with elements of vacuum, thus producing a high 
degeneracy. The low-degeneracy at the $^4 — $5 bound- 
ary is interpreted, from this vantage point, by the fact 
that hosts 2 do not accommodate hybrids 5 alone. 

By analogous reasoning based on the transcription ^ 
we can cast Eqs. | |11[) fo r the particles excited from the 
plateau state (Sec. II D ) in the form ( 13 ) with the right- 
hand sides changed into 



-2Kj Kj-K D -H Kj-Kd+H 



(21) 



and the associated grand partition function in the form 

JV 



Z= (l + w 3 )e- K '- K ^ 2 - 3H / 2 



(22) 



The calculation of the (N m ) produces the expected re- 
sults. The region with a low density of particles of any 
kind at low T is now $ p + . Regions $f+ , ^>s have high 
densities of compacts m = 1,2, respectively, region <E>p_ 
has a high density of tags to = 3, and regions §a,&p- 
have high densities of hosts m = 4, 5, respectively. In this 
case the low-degeneracy of the Q>a — 3>s phase boundary 
is explained by the observation that the negative-energy 
compacts to = 2 and the negative-energy hosts to = 4 do 
not interlink directly, only via (zero-energy) elements of 
pseudo- vacuum. 

From the expressions for Z and the (N m ) (all in four 
versions) a great many statistical mechanical results can 
be derived. However, the focus of this paper is on the 
the identification of functions associated with particles 
of given structures in specific environments. 



Equations (11 1 for the particles excited from the spin- 



polarized state[Sec. II C I can be simplified into Eqs. ( 17) 
with the right-hand sides changed into 



Kj+Kd+H 



-2H 



1. 



e -2Kj 
„2K, 



(19) 



The close relationship originates in the transcription pi) 
between the motifs of the two sets of particles. The grand 



IV. ENTROPY OF MIXING 

Depending on whether or not the pseudo-vacuum of 
choice is a state of lowest energy, the Ising chain ([I]) inside 
any of the six regions in the parameter space of Fig. [ija) 
can be interpreted (at low T) either as a dilute gas of sev- 
eral species of (positive-energy) particles or as a liquid of 
predominantly one species of (negative-energy) particles 
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with low concentrations of higher-energy particles added. 
The physical processes that take place when the Hamil- 
tonian parameters are gradually changed across one or 
the other phase boundary allow for multiple interpreta- 
tions in terms of particles from the four sets {ttia}, {ms}, 
{mp + }, {mp + }, introduced in Sec.|ll]and elements of the 
associated pseudo-vacua. 

We now consider the simplified phase diagram pertain- 
ing to ft = as shown in Fig. [TJb) . The singlet phase <&s 
and the antiferromagnetic phase §a are the pseudo- vacua 
of particles from the sets {ms} and {m^}, respectively. 
The ferromagnetic phase §f represents the degenerate 
phases $f±- One ground-state vector is the pseudo- 
vacuum for particles from set {m F+ } and the other is a 
solid of zero-energy particles from that set. The plateau 
phases $p± are not realized at h = 0. 

Along the two phase boundaries at D > 0, the degen- 
eracy of the ground state is low (threefold) and along the 
third phase boundary (at D < 0) it is exponentially high. 
This qualitative difference has a natural explanation ex- 
pressible in terms of particle functions. 



A. Surfactants 

We begin by examining the physics near the <!>a ~ &s 
boundary simultaneously from two vantage points. Phase 
<&A is the pseudo- vacuum of particles from set {m^} or a 
solid of hybrids mj = 7 (tlt> ItD- Conversely, phase <&s 
is the pseudo- vacuum of particles from set {ms} or a solid 
of tags ttia = 3 (00). At the phase boundary, the lowest 
excitations are separated by AE = J from the threefold 
ground state and contain exactly one host particle act- 
ing as a pair of barriers between one segment of tags or 
hybrids on the inside and one segment of pseudo- vacuum 
on the outside. Viewed from one vantage point, the host 
is among m A = 4, 5, 6 (t ° fj i I- 1 I, i t) and 
accommodates tags mA = 3. The other vantage point 
has one host ms = 1,2,3,4 (o f o,o 1 0,0 fj, 0,0 J/J- o) 
accommodating hybrids ms = 7. 

When we switch vantage points the host particle is 
turned inside out and the tags or hybrids inside become 
elements of pseudo-vacuum outside. One product state 
of this kind is 

1 : -n t°°---°°+ n-, > 

vac. host+tags vac. 

= I - ^W- - -tlt° °f - > (23) 

vac - host+hybrids vac - 

with host mA = 6 and tags mA — 3 embedded in the 
antiferromagnetic pseudo-vacuum or host ms — 4 and 
hybrids mj = 7 embedded in the singlet pseudo-vacuum. 

From either vantage point the host assumes the func- 
tion of a surfactant between elements of vacuum on one 
side and tags or hybrids on the other side. The absence 
of direct links between tags/hybrids and elements of vac- 



uum implies an infinite surface tension. The presence of 
surfactants reduces it to a finite value. 

At the phase boundary, the energy densities inside 
and outside the host are equal. Moving away from the 
phase boundary raises the energy density of one seg- 
ment and lowers that of the other, making the growth 
of one segment at the expense of the other energetically 
favorable. This process involves either the creation of 
negative-energy particles or the annihilation of positive- 
energy particles from the tag or hybrid categories. 

To get the transition started, one host particle (of posi- 
tive energy) must first be created. Its energy at the phase 
boundary is equal to J and decreases on a path into the 
phase with opposite pseudo- vacuum. Eventually it goes 
negative, producing a hysteresis effect. 

What happens when the phase boundary D = J is 
crossed at low but nonzero T evokes a metamorphosis 
reminiscent of an M. C. Eschcr print. As the tags = 3 
soften energetically, their population explodes. The hosts 
of the soft tags expand in the process. The number of 
thermally excited hosts remains low. The segments of 
pseudo- vacuum between the hosts shrink to sizes compa- 
rable to the expanded hosts. 

A change of perception identifies the close-packed soft 
tags of set {mA} as segments of pseudo- vacuum for set 
{ms} and segments of pseudo- vacuum for set {to a} as 
close-packed soft hybrids of set {ms}- The hosts from 
set {mi} break apart and reassemble into hosts from set 
{ms}- Each left half of an old host pairs up with the 
right half of the closest old host to its left to form a new 
host. 



B. Immiscible liquids 

Now we view the physics near the same $^4— $5 bound- 
ary from a third vantage point. We describe it in terms 
of particles from the set {mp+} whose pseudo- vacuum is 
not a lowest-energy state. Each of the three vectors that 
make up the physical vacuum at the phase boundary is a 
solid of one species of negative-energy particles: two vec- 
tors of hosts m F + = 2 (fit) an d one vector of hybrids 
m F+ = 5 (00). The energy per bond is — 2 J for both 
species. Hosts mp + = 2 do not accommodate hybrids 
mp + = 5. 

The lowest excited states contain one host mp + = 3, 4 
(t o |t,tl t) with energy per bond -5 J/3 at the 
phase boundary. These hosts interlink externally with 
hosts mp+ — 2 and internally with hybrids mp + = 5. 
They assume the function of surfactant by reducing the 
interfacial tension from infinity to a finite value. Vary- 
ing the parameters J, D across the phase boundary tips 
the energy balance from favoring hosts mp + — 2 to hy- 
brids mp + = 5 or vice versa. Any transition process 
that involves local spin flips cannot take place without 
the presence of a host mp + = 3, 4, which introduces the 
hysteresis effects noted previously. 

At low, nonzero T the thermodynamic state near the 
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FIG. 4: Average numbers (per site) (JVb), (N5) from the set 
{m,F+} and entropy (per site) S/ks versus J with D = 4 — J 
at fcsT = 0.25 (solid lines) and fcsT = 0.5 (dashed lines) in 
arbitrary energy units. The scale on the horizontal axis is for 
the lower T. The curves at the higher T are for 1.0 < J < 3.0. 



FIG. 5: Average numbers (per site) (iVi) + (iVa), (N3) from 
the set {tua} and entropy (per site) S/ka versus J with D = 
4 + J at k B T = 0.3 (solid lines) and k B T = 0.6 (dashed 
lines) in arbitrary energy units. The scale on the horizontal 
axis is for the lower T. The curves at the higher T are for 
-3.0 < J < -1.0. 



<f>A ~ ^>s phase boundary exhibits attributes of two im- 
miscible liquids. There are two species of particles with 
high population densities, the hosts mp + = 2 and the 
hybrids mp+ — 5. They coexist in an emulsion of sorts 
with higher-energy hosts mp + = 3,4 of much lower den- 
sity acting as emulsifiers. 

In Fig. [4] we have plotted (at two fixed values of T) 
the population densities (N2), (N5) of the particles that 
represent the two liquids versus J along a line that crosses 
the — &s boundary in a perpendicular direction. Also 
shown is the entropy per site as derived from ( |20[ ). 

Away from the phase boundary we have a fairly pure 
liquid of cither hosts m,F+ = 2 or tags mp + = 5. The 
different number densities are caused by the different par- 
ticle sizes. Particles of both liquids coexist in comparable 
numbers near the phase boundary (J = 2). 

The entropy has a maximum at the phase boundary. 
The key observation is that as T is lowered the coexist- 
ing populations of particles representing the two liquids 
remain high but the entropy is much reduced. It can be 
interpreted as the entropy of mixing of two immiscible 
liquids [32] . 

The physics near the <5>s~ &f boundary described from 
the vantage points of the four sets of particles is a vari- 
ation of the same theme. One significant difference is 
that the thermodynamic state as viewed from the vantage 
point of the particle set {toa} involves three immiscible 
liquids rather than just two. 

The particle species with high population densities at 
low T are the compacts ttia = 1,2 (tt,44-) at J < — D 
and the tags itia = 3 (00) at J > —D as shown in 
Fig. [5] Throughout region <!> p the negative-energy com- 
pacts remain segregated by higher-energy hosts vtia — 6 
(t 4-, 4- t) or elements of pseudo- vacuum (f4>, 4-t)- The 
entropy is very low. At the phase boundary the compacts 
coexist with tags tua — 3, but all three species remain 
segregated. The tags exist inside hosts and the compacts 



outside. 

In addition to hosts ttia = 6, hosts ttia = 4, 5 
(j" o i,], o i) are present in comparable numbers. The 
particle functions that segregate compacts from each 
other are different from those that segregate compacts 
from tags. The shallow peak in the entropy curve at low 
T is a dim signal of the increase in emulsion components 
near the phase boundary. 



C. Miscible liquids 

The physics near the $,4 — phase boundary is qual- 
itatively different. From the vantage points of the adja- 
cent phases the two vectors of phase &a are the pseudo- 
vacuum of {itia}, composed of elements tl-il-ti or two 
solids composed of hosts mp + = 2 (t4-t)> whereas the two 
vectors of phase &p are the pseudo- vacuum of {rrip + }, 
composed of elements ft, plus a solid of tags mp + = 6 
(44-) or two solids, one composed of compacts rriA = 1 
(ft) and the other of compacts = 2 (H). 

The low degeneracy of the pure phases is explained by 
the fact that the particles and elements of pseudo- vacuum 
involved in each pair of vectors do not mix. However, 
when both pairs of vectors become degenerate (at the 
phase boundary) the options for particles or elements of 
pseudo-vacuum to interlink mutually increase dramati- 
cally, producing a 2 Ar -fold degeneracy, which comprises 
all vectors \<J\---<jn) with oi =t,4-- Any such state 
can be interpreted as either a combination of compacts 
niA — 1,2 (1i\44-) and elements of associated pseudo- 
vacuum (t4->4-t) or as a combination of hosts mp + =2 
(tit)? tags mp + = 6 (4-4-), and elements of associated 
pseudo- vacuum (ff). 

From the vantage point of particle set {777,5}, the physi- 
cal vacuum at the $a~ &f phase boundary consists close- 
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packed configurations of particles with the lowest nega- 
tive energy per bond. They are the hybrids mj = 5, 6, 7 
(tt;-W;t-jT,-i-tl)- Their ability to interlink produces all 
2 N vectors of the ground state at the phase boundary. 
Varying the parameter J across the phase boundary at 
fixed D < enables the transition to proceed via local 
spin-flip processes without energy barrier. There are no 
hysteresis effects in this case. 

At low, nonzero T the thermodynamic state near the 
phase boundary can be interpreted as that of three liq- 
uids. The species of particles with high population den- 
sities are the hybrids ms = 5, 6 at J < and the hybrids 
m s = 7 at J > 0. The state in phase region $ F is akin 
to that of two immiscible liquids as already described (in 
Sec. IV B I from the vantage point of set {wia}. However, 
here the two relevant species are hybrids. In region $^4, 
by contrast, we have a one-component liquid, composed 
of hybrids ms = 7. The crucial difference in this scenario 
compared to the one described near the $s — <&f phase 
boundary is that the hybrids mg — 7 not only interlink 
with both species of hybrids, ms — 5, 6, but they also 
mediate their mixing by acting as a surfactant with zero 
intcrfacial energy. 

Figure [6] shows the population densities of hybrids 
ms — 5, 6, 7 near the &a — ®f phase boundary and also 
the entropy, all at two different temperatures. We ob- 
serve a strikingly different behavior of the entropy curve 
in comparison to its behavior across the other two phase 
boundaries (shown in Figs. [4j |5j. Here the particles 
with dominant populations behave as in miscible liquids. 
Again there is a region where these populations coexist 
in comparable high densities at both the higher and the 
lower temperature. Unlike in the previous two scenarios 
the entropy stays high at the lower temperature because 
the particle populations are not separated by barriers in 
the form of higher-energy host particles acting as surfac- 
tant with significant interfacial energies. 
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FIG. 6: Average numbers (per site) {N5) = {Ne}, {N7) from 
the set {ms} and entropy (per site) S/ks versus J at D — —2 
at ksT = 0.15 (solid lines) and ksT = 0.5 (dashed lines) in 
arbitrary energy units. The scale on the horizontal axis is the 
same for all curves. 



V. STRUCTURE AND FUNCTION 

In this project we started out with interacting Ising 
spins on a periodic chain. The model Hamiltonian ([TJ 
causes the self-assembly of particles with various struc- 
tures that (i) remain free of interaction energies even at 
high densities and (ii) are robust against changes in the 
Hamiltonian parameters. 

We have identified four sets of particle species excited 
from four different ordered Ising states. From there on 
we are free to remove the scaffold of the Ising model and 
assign arbitrary energies to all particles species |28j . 

Each particle species from a given set has a distinct 
structure. That structure is encoded in one or several mo- 
tifs. The structure of a particle is only meaningful in re- 
lation to the structure of the associated pseudo-vacuum. 
Particles tua — 2, ms — 6, and mp + — 3, for exam- 
ple, all have the same motif, H, but represent different 
structures. They belong to different species, excited from 
different pseudo-vacua. The statistical interaction speci- 
fications of all three species are different. 

We have classified all particle species into four cat- 
egories according to structural features. These features 
represent specific particle functions. Compacts have only 
external links. They interlink with elements of pseudo- 
vacuum and with external links of compacts and hosts. 
Hosts have external and internal links. They interlink ex- 
ternally with elements of pseudo-vacuum and with com- 
pacts and hosts. They interlink internally with tags or 
hybrids. Tags and hybrids interlink with internal links 
of hosts and, selectively, with each other. These limited 
options for particles to interlink mutually may be inter- 
preted as a manifestation of stereospecific binding with 
zero binding energy. 

The capacities of hosting and being hosted are func- 
tions that some particle species have and other species 
lack. Some functions of particle species depend on the 
circumstances as is the case in two examples encoun- 
tered in Sec. |IV| (i) In a situation where two particle 
species that do not interlink directly have equal negative 
energy density, lower than the energy density of any other 
species, then they form two immiscible liquids a low T. 
Segments of one liquid are separated from segments of 
the other. The function of surfactant is assumed by a 
third particle species or by elements of pseudo-vacuum 
with higher density, (ii) If the circumstances change in 
such a way that the energy density of the surfactant par- 
ticle species becomes equal to that of the particle species 
representing the two immiscible liquids then the surfac- 
tant particle changes its function into a miscibility agent 
of sorts. 

Modifications in the model Hamiltonian affect the 
nature of particles. We have already seen that variations 
in the Hamiltonian parameters J, D, h cause changes in 
some functions of particles even though their structures 
remain invariant. 

Additional commuting terms in the Ising Hamiltonian 
such as a next-nearest-neighbor coupling or a bond al- 
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ternation make it necessary to find new sets of particles 
with longer motifs that interlink differently and exhibit 
different statistical interactions. Ising Hamiltonians with 
spin s > 1 pose similar challenges apart from the need 
for motifs with new symbols |33) . 

Non-commuting terms in the Hamiltonian are apt to 
affect the particles introduced here more drastically and 
in multiple ways. The particles are likely to move along 



the lattice, scatter off each other, and decay into other 
particles. Under special circumstances as mentioned in 
Sec. [TJ the particles do not decay and experience only 
elastic two-body scattering. Such particles can be char- 
acterized by motifs in momentum space [3] . A taxonomy 
of particles according to structures and functions in this 
wider arena is, of course, a more challenging project. 
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